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ABSTRACT 

We present observations of two strongly lensed z ~ 5 Lyman-a Emitting (LAE) galaxies that 
were discovered in the Sloan Giant Arcs Survey (SGAS). We identify the two sources as SGAS 
J091541+382655, at z = 5.200, and SGAS J134331+415455 at z = 4.994. We measure their AB 
magnitudes at {i, z) = (23.34 ± 0.09, 23.29 ± 0.13) mags and (i, z) = (23.78 ± 0.18, 24.241°;^^) mags, 
and the rest-frame equivalent widths of the Lyman-a emission at 25.3 ± 4.lA and 135.6 ± 20.3A for 
SGAS J091541-I-382655 and SGAS J13433H-415455, respectively. Each source is strongly lensed by 
a massive galaxy cluster in the foreground, and the magnifications due to gravitational lensing are 
recovered from strong lens modeling of the foreground lensing potentials. We use the magnification to 
calculate the intrinsic, unlensed Lyman-a and UV continuum luminosities for both sources, as well as 
the implied star formation rates (SFR). We find SGAS J09154H-382655 and SGAS J134341-K415455 
to be galaxies with {hhy-a, L(7v) < (0.6L|^j^_„, 2L^y) and {IjLy-a, Lc/y) = (0-5L^j^_^, 0.9L^y), 
respectively. Comparison of the spectral energy distributions (SEDs) of both sources against stellar 
population models produces estimates of the mass in young stars in each galaxy: we report an upper 
hmit of Mstars < 7.9l|j X 10^ Mq/iq j for SGAS J091531-I-382655, and a range of viable masses for 
SGAS J134331+415455 of 2 x 10** Mq/iq;^ < M^tars < 6 x 10^ Mq/iq 
Subject headings: gravitational lcnsing:strong ~ galaxies: high-redshift galaxies 



1. INTRODUCTION 

Understanding the evolution of galaxies ~ especially 
the first generation of galaxies - remains one of the 
most important topics in astrophysics and cosmology. 
Populations of galaxies in the distant universe are 
identified over a wide range of wavelengths, including 
Distant Red Galaxies (DRGs), Ultra Luminous Infra- 
Red Galaxies (ULIRGs) and Sub-Millimieter Galax- 
ies (SMGs). Many efforts to study the properties of 
high redshift galaxies at optical wavelengths focus on 
two distinct classes selected by their rest-frame UV 
properties: 1) Lyman-Break Galaxies (LBGs) and 2) 
Lyman-a Emitters (LAEs). LBGs are selected via deep 
wide-band photometry, identified by the 'Lyman limit' 
continu um break that app ea rs at 9 12 A in the r est 
frame (jSteidel et al.l Il996allbl : iLowenthal et al.l 119971 ) - 
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though for sources at higher redshift this spectral break 
moves redward, approaching 1 216A in the rest frame 
due to the Lyman-a forest (jSteidel fc Sargent 119871 : 
lRauch|[T998D absorption by intergalactic neutral hydro- 
gen. LAEs ar e selected by either nar row-band i maging 
(iCowie fc H"i]|T998: R hoad s et al.ll200'0. 2(303: Aiiki etal] 
120041: IGawiser et al.1^^0061: lYamada et al.lT2005l) o r blind 
spectroscopy (|Kurk et all 120041: ISawicki et al.l I2008D 
tuned to detect Lyman-a line emission redshifted into 
near-ultraviolet, optical, or near-infrared wavelengths. 
Over the past decade large samples of LBGs and LAEs 
have driven studies of star-forming galaxies at z > 2.5. 

Surveys for LBGs and LAEs are efficient for collect- 
ing statistical samples of high-redshift galaxies, but at 
z > 3 they produce objects which are generally too faint 
to have their galactic continuum emission studied spec- 
troscopically. The standard approach for studying the 
properties of these galaxy samples relies on stacking the 
photometric signal from many objects and fitting the ob- 
served mean SED against a variety of stellar population 
synthesis models in order to constrain parameters such as 
the ages and masses of the underlying stellar popul ations. 



as well as the amo unt of dust extinction (Shaple v et al 
2003; Ch ary et al.l [ 2005: Pirzkal ct al. 2007; iLai et al 



| 2007i: iFinkelstein et al...200& . ,2009a; Nilsson et al.ll2009l : 
lYabe et al.ll2009( )" In principal, stellar population syn- 
thesis modeling can also provide information about dust 
properties (i.e. the shape of the dust law) and mctal- 
licity, but even the stacked SED signal at z > 3 is in- 
sufficient to constrain these additional parameters with 
much confidence. Broadly speaking, galaxies selected as 
LBGs are believed to sample more massive star-forming 
galaxies with an underlying older stellar population, and 
possibly higher dust content, while LAE selected galaxies 
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tend to be lower mass galaxies wi t h low metallicities and 
very little dust (lGiavaliscol l2002t iVenemans et alj 120051: 
iGawiser et al.|[2007D . Hubble Space Telescope imaging 
studies of high redshift LAE galaxies imply that these 
sources are compac t, and likely either disk-like or ir- 
regul ar in struct ure (iPirzkal et alJ l2007l: iTaniguchi et_a]J 
120091 ). Recently iFinkelstein et"ar ( 2009^ modeled indi- 
vidual SEDs of 14 bright z ~ 4.5 LAEs from the Chandra 
Deep Field South and found a broad range in stellar pop- 
ulation age, stellar mass, and dust extinction, which sug- 
gests that stacking SED analyses of high redshift galaxies 
may not be the best approach. 

The main hurdle involved in studying any high redshift 
source is the general lack of signal. Distant galaxies are 
faint and therefore difficult to detect, and those which 
are identified are rarely - if ever - amenable to detailed 
follow-up. Furthermore, those sources which are bright 
enough to be studied individually are drawn from the 
extreme bright tail of the luminosity function of high 
redshift galaxies, and are therefore not necessarily repre- 
sentative of the bulk of the populations. In this paper we 
present two serendipitously discovered, strongly lensed 
high redshift galaxies: SGAS J091541-I-382655, spectro- 
scopically confirmed at z = 5.200 ± 0.001, with r^B — 
24.68 ± 0.25, iAB = 22.92 ± 0.09 and zab = 22.75 ± 0.13 
mags, and SGAS J134330-I-415455 spectroscopically con- 
firmed at z = 4.994 ± 0.001, with tab > 25.47 mags, 
iAB = 23.36 ± 0.18 mags and zab = 23.70^0 ijg mags. 
Both objects have riz colors and magnitudes that sat- 
isfy selection criteria for r— band dropouts in z ^ 5 
dropout surveys, as well as Lyman-a equvialent widths 
(see Section 2) sufficiently large to be selected in sur- 
veys for Lyman-a excess. At z > 5 these objects are 
the two brightest LAEs in the literature to date, and 
both sources are projected on the sky within 30" of the 
cores of confirmed strong lensing galaxy clusters. This 
means that the sources - once corrected for the lensing 
magnification - are intrinsically much fainter than the 
observed flux implies, and therefore provide a rare oppor- 
tunity to study individual LAE properties at the fainter 
end of the luminosity function. There is a small but 
growing number of magnified galaxies at high redshift 
that are excellen t candidates for high resolution spectro - 
scopic follow-up (|Koester et al.|[2010l: Imivts et al.|[20Tol ): 
some of these galaxies have been observed in detail at 
optical and n ear-infrared wavelengths, including cB58 
(IPettin i et al.i r200a). "the 8 O 'clock Arc" ((Allani et al.< 
I2007t I Finkelst ein et all I2009bf), and "the Cosmic Eye " 
(jSmail et al.u2007l: ISiana et al.ll2009l: iQuider et al1[2010l ) . 



The two galaxies discussed in this paper are the first 
z ^ 5 galaxies that present similar opportunities for de- 
tailed study via follow-up spectroscopy. 

Where necessary we calculate cosmological distances 
assuming a flat cosmology with Hq = 70 km s~^ Mpc~^, 
and matter density f2j\/ — 0.3. All magnitudes are AB. 

2. OBSERVATIONS 
2.1. Data 

The two sources presented here were first identified as 
r-band dropouts in gri imaging of two different strong 
lensing galaxy clusters, and subsequently confirmed by 
spectroscopy to have strong Lyman-a emission features. 
The two galaxy clusters were identified as part of the 




Fig. 1. — CMOS griz (from left to right) 8" X 8" cutout images 
centered on SGAS J091541-I-386254 and SGAS J134331-I-415455. Both 
objects exhibit obvious drop-out behavior from the i— to r — bands. 
The 2: — band detection of SGAS J091541+386254 appears at lower sig- 
nificance than its counterpart for SGAS J134331+415455 because its 
z data were taken during bright time and the sky background noise is 
~ 30% higher. 



Sloan Giant Arcs Survey (SGAS; iHennawi et"ani2008| ). 
a blind survey for strong lensing systems in optically se- 
lected massive clusters at 0.1 < z < 0.6 detected via 
the R ed-Sequence Cluster algorithm ()Gladders &: Yee I 
I2OOOD a dapted to run on the Sloan Digital Sky Survey 
(SDSS; lYork et all I2000D public data release catalogs. 
Strong lensing clusters are identified by visual inspec- 
tion of imaging in 5— band on 2m to 4m-class telescopes 
for 600s in < 1" seeing, and the most spectacular sys- 
tems have been followed-up with multi-band imaging and 
spectroscopy on 8m-class telescopes. One of the clusters 
discussed here, SDSS J1343+4155, appears in a recent 
small sam ple of stron g lensing clusters discovered in the 
SDSS bv Diehl et all ([2009). The other cluster, SDSS 
J0915+3826, does not appear in any prior published 
work. Imaging and spectroscopic observations were con- 
ducted with the Frederick C. Gillett Telescope (Gemini 
North) between the months of 2008 February and 2008 
July. The CMOS imaging observations were pre-imaging 
conducted for the purpose of spectroscopic mask design. 
Both imaging and spectroscopy were part of Gemini pro- 
gram GN-2008A-Q-25. The primary goal of the spectro- 
scopic observations was to obtain redshifts of arcs to fa- 
cilitate strong lensing modeling. We briefly summarize 
the spectroscopy here, and refer the reader to Bayliss et 
al. (2010), in preparation, for additional details. 

Gemini North/CMOS gri photometry for both cluster 
fields are derived from 2xl50s dithered exposures, which 
were executed in queue mode on February and March 
of 2008. These pre-imaging data were used for spec- 
troscopic mask design and target prioritization. CMOS 
z-band observations consisted of 6x1 80s dithered expo- 
sures which were scheduled as follow-up, primarily in or- 
der to measure the continuum fiux of the LAEs, and were 
executed in queue mode in February of 2010. The CMOS 
images were reduced using the Gemini IRAlfl package. 

In addition to the GMOS-N photometry and spec- 
troscopy, we obtained near-infrared (NIR) imaging of the 
two lensing clusters in the z'JH filters with the Near- 
Infrared Camera and Fabry-Perot Spectrometer (NIC- 
FPS) of the 3.5m telescope at the Apache Point Obser- 

IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 
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vatory (APO) in New Mexico. The detector is a Rockwell 
Hawaii 1-RG 1024x1024 HgCdTe device with a pixel scale 
of 0.273" pixel"^ and a 4.58 x 4.58 arcmin unvignetted 
field of view. The APO/NIC-FPS z' data differs from 
the Gemini/GMOS z significantly due to the different 
wavelength responses of the NIC-FPS and GMOS de- 
tectors. The effective wavelengths of the two filters are 
offset by ~ lOOOA, and we use the prime (') through- 
out this paper to distinguish the APO/NIC-FPS z' from 
the Gemini/GMOS z. The NIR data were taken on 3 
different nights in the Winter and Spring of 2009. The 
conditions during the observing nights varied, with sub- 
arcsecond seeing for the SDSS J1343-I-4155 z'-band and 
both J-band images. The SDSS J0915-h3826 z'-band and 
both H-band images were taken in ^ 2" seeing. The ob- 
servations consisted of 5-point dithers around a 40" box 
and were reduced, registered and stacked using a custom 
IRAF pipeline. Total exposure times are 5400s, 6960s 
and 3500s for 2'-, J-, and /f-band observations of SDSS 
J0915-I-3826 and 7800s, 5160s and 2625s for z'-, J-, and 
77-band observations of SDSS J1343-I-4155. 

The SDSS was used to calibrate the four optical bands 
and the Two Micron AU Sky Survey (2MASS) was used 
to calibrate the NIR observations. Prior to making 
any photometric measurements we first transform im- 
ages in all bands to the reference frame of the z-band 
image, 0.1454" pixel" ^ (GMOS-North detector, binned 
2x2). We then construct an empirical, normalized point 
spread function (PSF) for each image based on a well- 
defined, non-saturated reference star. We create photo- 
metric apertures by drawing a ridge line that covers the 
high-redshift LAE and convolving it with the appropri- 
ate PSF for each image. Apertures are defined by an 
equivalent radius, which corresponds to the radius of a 
circular aperture that goes out to the same isophot. We 
make the final magnitude measurement using a detailed 
sequence of sky subtraction and outlier masking steps; 
first we subtract a general sky measurement and then 
compute the median pixel value and standard deviation 
inside annuli of fixed width at increasing radial distances 
from the source. Far enough out, these median values 
converge to zero for an accurate sky subtraction. We 
average the median values at large radii and subtract 
this average from the image to correct the general sky 
subtraction. Outliers are defined as pixel values that de- 
viate by more than 5cr from the median in the respective 
annuli, and are replaced by the median value plus an ap- 
propriate noise term. The final magnitude is measured 
at an equivalent radius of twice the FWHM of the PSF 
and corrected to an equivalent radius of 6" based on the 
curve of growth of the PSF. By defining the radius as 
twice the FWHM, we make sure our apertures always 
cover the same physical region on the sky. In the case 
of SGAS J091541-I-382655, a foreground galaxy lies very 
close to and partiall y on top of the L AE galaxy. We use 
the GALFIT package (|Peng et al.ll2002[ ) to fit a Sersic pro- 
file to this galaxy in the imaging data where the LAE has 
no measureable flux, and then scale the galaxy model by 
the peak flux in each band and subtract it out before 
measuring the LAE magnitudes. 

When the LAEs are not detected (in gJH for SGAS 
J091541+382655and grJH for SGAS J134331+415455), 
we measure a limiting magnitude, derived from the total 



sky noise in the aperture. We define a sky noise per 
pixel (Tsfcy as the standard deviation of our overall sky 
measurement plus a contribution from the poison error 
made in this measurement. A pixel value of 2(Jsky is 
added in quadrature for each pixel in our final aperture 
with an equivalent radius equal to twice the FWHM. 
The limiting magnitudes are also aperture corrected to 
an equivalent radius of 6". This method constrains the 
source to be fainter than the limiting magnitudes at 95% 
confidence in the relevant bandpasses. 

All spectroscopic observations were carr ied out us- 
ing t he Gemini Multi-Object Spectrograph (jHook et al.l 
120041) using custom slitmasks that were designed to tar- 
get lensed sources based on their color, morphology, and 
positoin in the GMOS imaging. After targeting all of 
the arc candidates, any remaining slits were placed on 
cluster members, easily identified by their red sequence 
colors. Spectra were taken using the macroscopic nod- 
and-shuffle (N&S) mode available on GMOS. The reasons 
for using N&S are threefold. First, many of the emission 
or absorption features that are used to determine galaxy 
redshifts in the range z = 1.0 — 3.0 characteristic of the 
giant arcs, are in the redder part of the optical where 
sky lines are problemati c, and NfcS facilitates more ac- 
curat e sky-subtraction ([Glazebrook fc Bland-HawthornI 
l2001f ). particularly at low spectral resolution. Second, 
N&S sky-subtraction allows us to use very small 1" x 1" 
microslits that can be densely packed into the cluster 
core 30"), allowing us to targe t as many arcs, arcle ts, 
and cluster galaxies as possible (jGilbank et al.l l2008f ) . 
Third, as we are primarily interested in a number of ob- 
jects around the cluster center, the density of objects 
and the limited field of interest are perfect for block- 
shuffling N&S observations. Spectra were taken with the 
R150_G5306 grating in first order which gives a disper- 
sion of 3.5A per pixel (binned spectrally by two), with 
six pixels per resolution element resulting in a spectral 
FWHM ~ 940 km s'^ Although the R150 grating of- 
fers broad spectral range from the atmospheric cutoff to 
A > 1/im, the drop in sensitivity at the blue and red 
extremes, due both to the GMOS CCD and the R150 
grating efficiency, results in effective spectral coverage of 

- 4000 - 95OOA . 

The N&S technique employed in our program is non- 
standard in that it involves a nod distance on the sky 
that is half the size of the macroscopic shuffie. The mask 
is designed to incorporate two submasks, each of which 
is a set of slits covering an area one third the size of the 
detector. Slits for the two submasks overlap on the sky in 
an area that is one sixth the size of the detector (because 
the nod distance is set to half the shuffie distance). This 
design allows us to place science slits for the primary 
target - the core of a strong lensing cluster in this case 

- on both submasks and obtain useful spectra for the 
entire duration of the N&S exposure, whereas standard 
macroscopic N&S results in science spectra for only half 
of the total exposure time. Our integration times were 
2400s resulting in a 1200s effective integration for each 
of the two submasks. Two exposures were taken for each 
target. Thus if an arc was targeted on both submasks 
(typical for the most prominent arcs) the total integra- 
tion time was 4800s. N&S facilitates straightforward sky 
subtraction by differencing two sections of the detector 
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- each 1/3 the size of the full detector. The spectra 
presented here were wavelength calibrated, stacked, ex- 
tracted and analyzed using a custom pipeline based on 
the XIDL software packageQ. 

2.2. SGAS J091541+382655 

SDSS J0915+3826 is a new strong lensing cluster that 
was discovered in the SGAS - discussed above. We 
measure redshifts for 16 cluster members in our Gem- 
ini/ GMOS spectroscopy and combine this new data with 
a redshift for the BCG measured in the SDSS to identify 
a mean cluster redshift oi z = 0.397 and a velocity dis- 
persion of 870 ± 169 km . The most prominent strong 
lensing feature around this cluster is a bright blue arc 
that we identify as a background galaxy at z = 1.501 
from 0[II]A3727A, CII]A2328A, CIII]A1907, 1909A in 
emission, as well as AIIIA1855, 1863A in absorption. Pre- 
imaging of this cluster reveals a source near the cluster 
core that exhibits a dramatic drop in flux from the i- to 
r-band, as shown in Figure [TJ With Gemini/GMOS-N, 
in griz, we measure this source at gAB > 26.07, tab = 
24.86±0.25, iAB = 23.34±0.09, and zab = 23.29±0.13. 
Slits were placed on this source in both submasks for our 
N&S observation of this cluster and the spectra exhibit 
bright emission at 7539A which we interpret as Lyman-a 
A1216A at 2 = 5.200. 

The LAE galaxy in this case is located on the sky 
very near to a foreground galaxy, and so we took care 
to account for possible contamination in the LAE spec- 
trum by light from the galaxy that falls into the slit, as 
this could mimic a continuum break in the source spec- 
trum. The separation between the LAE box slit and 
the potentially contaminating foreground galaxy is 2.5". 
Examination of the radial profile of the potential con- 
taminant in the pre-imaging observation shows that the 
i-band fiux within the LAE slit due to the galaxy is on 
average 1.5x10"'^ of the peak flux of that galaxy. The 
mean i-band flux of the LAE within the spectroscopic 
aperture is 5.1x10"^ compared to the peak of the con- 
taminating galaxy, and thus assuming comparable see- 
ing between the imaging and spectroscopic observations 
suggests that this object's contribution to the i-band flux 
within the LAE aperture is approximately 3% of the total 
light. Note that given the separation between the con- 
taminant and the slit the seeing of the spectroscopic and 
imaging observations would have to be grossly different 
to produce significant contamination; such a mismatch 
is neither expected due to execution of the observations 
in specified conditions in queue mode, nor suggested by 
the spatial width of stellar spectra within the alignment 
boxes of the spectroscopic observations. 

The relevant part of the extracted spectrum is dis- 
played in Figure [5J and despite the low dispersion grating 
used for our observations the line is measurably asym- 
metric. We also detect continuum signal redward of the 
emission line and no significant signal blueward of it. 
This source has a strong photometric r-i break - even 
after subtracting the emission line flux from the z-band 
photometry. It is also very blue in i-z, and is not de- 
tected in NIC-FIPS imaging to a 2-a limiting magnitude 
of Jab > 22.49, indicatingg that the continuum spec- 

^ http:/ /www. ucolick.org/~xavier/IDL/index. html 
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Fig. 2. — GMOS spectrum of the Lyman-CK emitter behind cluster 
lens SDSS J0915+3826. The LAE spectrum plotted over a range of 
~ 7150 — 8350A is smoothed to a scale comparable to a the spectral 
resolution of 940 kms^^ in order to emphasize the low-significance can- 
didate absorption features, with the composite Lyman Break Galaxy 
absorption+emission spectrum from Shaplcy ct al. ( 2003) overplottcd 
in green. Dashed vertical lines indicate the location of prominent UV 
continuum metal absorption lines at the LAE redshift, which coincide 
with the low-significance features in our continuum spectrum. The 
Lyman-a emission feature is also shown inset, without smoothing. A 
typical spectral flux error bar is displayed on the right side of the figure, 
in the whitespace above the continuum. 

tral slope redward of the emission feature is blue. The 
blue continuum redward of the emission and the strong 
r-i break - along with the non-detection in g, and the 
asymmetry of the line confirm the interpretation of the 
source as a Lyman-a emitting galaxy at high redshift. 
The most likely candidate emission lines which confuse 
searches for high-redshift LAEs are 0[II]A3727A, [OIII] 
A5007A and H-a, but for each of these other candidate 
lines we would expect accompanying emission features to 
fall within the wavelength coverage of our spectroscopy. 
We measure the equivalent width of the Lyman-a emis- 
sion by fitting a simple power-law to the continuum spec- 
trum redward of the emission line over a wavelength 
range of 7575 — 8275A, and subtracting the underlying 
continuum flux from the emission line. Over the rela- 
tively small wavelength range used the continuum level 
is consistent with a constant value. We note that our 
spectral flux calibration uses a standard from the Gem- 
ini archive to correct roughly for the detector sensitivity 
as a function of wavelength, but is not suitable for mea- 
suring absolute fluxes from the spectral data. Taking 
into account the underlying continuum and the source 
redshift we measure a rest-frame EW for the Lyman-a 
emission of 25 ± 4A, and combine the EW measurement 
with our i— band photometry to calculate an observed 
Lyman-a line flux of 2.1 ± 0.7xl0~^^ ergs s~^ cm~^. 

2.3. SGAS JI3433I+415455 

The second lensed LAE presented here is located near 
the strong lensing cluster SDSS J1343+4155, which was 
previousl y identified as a s tron g lensing clu ster in the 
SDSS bv lWen et al~l ^OM) and iDiehl et al.l pOOaV We 
combine our Gemini/GMOS spectroscopy with two clus- 
ter member redshifts from the SDSS to measure a mean 
cluster redshift oi z = 0.418 and a velocity dispersion 
of 1077 ± 266 km s^^ from 7 cluster member galax- 
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Fig. 3. — GMOS spectrum of the Lyman-a emitter behind cluster 
lens SDSS J1343+4155. The spectrum is presented identically to the 
spectrum shown in Figure [2] except that it covers an observer-frame 
wavelength range of ~ 6900 - SIOOA. 



ies. The bright a rc around this len s is spectroscopi- 
cally identified by iDiehl et al.l (|2009| ) as a background 
galaxy at z= 2.091. Pre- imaging of this cluster reveals 
a source near the cluster core that exhibits a dramatic 
drop in flux from the i- to r-band - shown in Figure [TJ 
We measure this source at gAB > 25.97, tab > 25.64, 
iAB = 23.78 ± 0.12, and zab = 24.24^°;^^. We placed 
slits on this source in each of the two submasks for our 
spectroscopic N&S observations of this cluster. The spec- 
tra corresponding to these slits exhibit a bright emission 
line at 7289A, which we interpret as Lyman-a A1216A at 
z = 4.994. 

The extracted spectrum around the emission feature 
is displayed in Figure |3l The emission feature is sig- 
nificantly asymmetric and we measure continuum emis- 
sion redward of the line, but no continuum blueward. 
The source is undetected in NIC-FIPS imaging down 
to a 2-a limiting magnitude of Jab > 22.22. Simi- 
larly to the case of SGAS J091541+382655, this source 
has a blue i-z color and our NIR photometry implies 
that it has a blue continuum spectral slope redward of 
the emission feature. The most likely candidate emis- 
sion line which could be misinterpreted as Lyman-a in 
this spectrum is 0[II] A3727A, but we rule this out as 
a realistic interpretation based on the absence of H-/3 
A4862A and 0[III] A4960, 5007A, both of which should 
fall just within the spectral range covered by our data 
if the source were a very red galaxy at lower redshift. 
Other common contaminants in searches for high-redshift 
LAEs include [OIII] A5G07A and H-a, but in each of 
these cases we would also expect to observe other ac- 
companying emission features given the wavelength cov- 
erage of our spectroscopy. We measure the equivalent 
width of the Lyman-a emission by fitting a simple power- 
law to the continuum spectrum redward of the emission 
line over a wavelength range of 7400 - 8275A, and sub- 
tracting the underlying continuum fiux from the inte- 
grated flux from the emission line. Similar to the contin- 
uum fit for SGAS J091541-)-382655, the continuum level 
for SGAS J134331+415455 is consistent with a constant 
value. Taking into account the continuum and the source 



redshift we measure a rest-frame EW for the Lyman-a 
emission of 136 ± 20 A and an combine the EW mea- 
surement with our i— band photometry to calculate an 
observed Lyman-a line fiux of 2.1 ± 0.5 x 10^^® ergs s^^ 



3. ANALYSIS 

3.1. Lens Models and Intrinsic Source Properties 

Using the measured equivalent widths and i magni 
tudes we also calculate lensed isotropic Lyman-a line 
luminosities of Lij,_Q, 



5.9 ± 1.1 X 10^3 g^g s-i/i2 7 



and LLy-a = 5.4 ± 0.9 x lO^^ g-i/ig^ for SGAS 
J091541+382655 and SGAS J134331+415455, respec- 
tively. These are not the true isotropic luminosities be- 
cause both sources are lensed by foreground galaxy clus- 
ters, and are therefore significantly magnified. To mea- 
sure the intrinsic luminosities of the LAEs, we estimate 
the magnification due to strong lensing by the interven- 
ing clusters. The mass models are constr ucted using 
the p ublically available software, LENSTOOL () Julio et al.l 
Hop), with Monte Carolo Markov Chain (MCMC) min- 
imization in the source plane, and are shown in Fig- 
ure m For SDSS J1343+4155, we compute a simple 
mass model using as constraints the giant blue arc at 
z = 2.09, the LAE a,t z = 4.994, paired with a counter 
image candidate that we have identified. The cluster and 
the brightest cluster galaxy (BCG) are represented by 
pseudo isotherma l elliptical mass distributions (PIEMD; 
lJullo et all 120071 ). We allow all the parameters of the 
cluster halo to vary. For the BCG, we follow the light 
distribution for the positional parameters, and vary the 
core and cut radii and the velocity dispersion. Our best- 
fit model reproduces the locations and orientations of 
the observed lensed images, and is consistent with the 
measured velocity dispersion to within the uncertainty. 
Based on this model, the magnification at the location 
of the LAE is m ~ 12, which we use to calculate the 
intrinsic isotropic Lyman-a line luminosity for SGAS 
J134331+415455 to be LLy-a = 4.5 ± 0.7 x 10*^ erg 
s~^/ig 7. Assuming Case B recombination and taking the 
prescription from Kennicutt (1998) we calculate the star- 
formation rate (SFR) for SGAS J134331+415455 to be 
SFRiy.a = 4.1 ± 0.6 M0 yr^i/io;^. Applying a mangifi- 
cation correction to the apparent magnitude in z— band 
we find the intrinsic z— band AB magnitude to be 2; = 27. 
We can also estimate the SFR from the UV continuum 
flux density at 1500 A according to equation (1) from 
iKennicuttl (|199a ). which rcsuhs in SFRuv = 18 ± 3 Mq 

The strong lensing model for SDSS J0915+3826 is con- 
strained by the positions of the components of the giant 
arc. We do not identify a clear counter image for the 
LAE, and therefore do not use it as a constraint. We 
represent the lens with a single PIEMD and allow all its 
parameters to vary, which produces a velocity dispersion 
estimate that is similar to the value we measure from 
spectroscopy of cluster members (see Section 2.2). The 
critical curves for our best-fit model are plotted in Fig- 
ure [5l Since the available constraints are internal to the 
radial projection of the LAE, the predicted LAE magni- 
fication is highly uncertain. In particular, if we compute 
a set of models with parameters in the range allowed by 
the uncertainties, the location of the critical curve for 
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Fig. 4. — The field of SDSS J1343+4155 with tangential and radial 
critical curves for the bright main arc overplotted. One image of the 
lensed background LAE, spcctroscopically confirmed at z — 4.994, and 
its putative counterimage predicted by the lensing model are identified 
by white triangles. The center of the dark matter halo and BCG mass 
components are indicated by a square and cross, respectively. Other 
sources in the field — background objects and cluster members — arc 
labeled with their spectroscopic redshifts. Color images are created 
from Gemini+GMOS-North gri 300s pre-imaging data. 




Fig. 5. — The field of SDSS J0915+3826 with the tangential critical 
curve for the bright main arc overplotted. Other sources with spec- 
troscopic redshifts are also identified, including the lensed background 
LAE at z — 5.200. Color images are created from Gemini+GMOS- 
North gri 300s pre-imaging data. 



the LAE varies significantly and there is practically no 
upper limit for the implied magnification. Since in all 
cases the magnification is higher than ~ 10, we adopt 
this number as our working order-of-magnitude estimate 
and determine an upper limit for the intrinsic isotropic 
Lyman-a line luminosity for SGAS J091541-h382655 of 
Lli/~q < 5.9 ± 1.1 X 10 ergs s~^ft.Q 7 and a correspond- 
ing SFRiy.a < 5.3 ± 1.0 M0 yr-i/1,^7. Correcting the 



z— band magnitude for the lensing mangifciation we re- 
cover an intrinsic magnitude of z > 25.75 mags, from 
which we measure SFR f/y < 51 ± 7 Mq yT~^hQ]. using 
equation (1) from Kennicuttl" (|1998l ). 

We note that SFRlj,-q and SFRfyy differ signifi- 
cantly for both sources. For SGAS J134331+415455 
SFRiy-a is lower than SFRf/y by a factor of 4.5, 
which is consistent with previous studies finding that 
the SFR calculated from Lyman-a emission is often sys- 
tematically lower than other SFR metrics - such as 
UV continuum fi ux density at 1500 A - by a factor 
of ~ 5 or more (|Tapken et all [2001 . The disa gree- 
ment between the two SFR estimates is more extreme 
for SGAS J091541-I-382655, with SFRuv ~ 9x larger 
than SFRiy-a- In a study of LAEs at z = 3.1 iden- 
tified in the Chandra Deep Field South, iGronwall et alj 
(pool find on average, SFR[/y/SFRiy_o, ^ 3 with a 
large scatter about the mean, including several objects 
with SFR(7y/SFRLy-Q as large as we measure for SGAS 
J091541+382655. Significant differences between these 
two SFR calibrations are understandable considering the 
potential systematic errors involved in measuring the 
SFR from both UV continuum and Lyman-a emission. 
For example, resonant scattering of Lyman-a photons 
off of neutral hydrogen can supress observed Lyman-a 
emission if there is even a very small amount of dust 
in the source galaxy, thus producing an underestimate 
of SFRiy_a. Additionally, Kcnnicutt (1998) identifies 
several caveats associated with the SFRf/y calibration, 
including assumptions about the IMF and the timescale 
and manner of star- form ation (e.g. consta nt vs. burst). 
Specifically, equation 1 in lKennicuttl (|1998[ ) recovers SFR 
from the rest-frame UV flux by assuming a galaxy has un- 
dergone continuous star formation for 10* years or longer, 
and originates from a Salpeter IMF with mass limits of 
0.1 and 10 M0. If the properties of the underlying stellar 
populations in either of these galaxies differ signifi cantly 
from t he parameters assumed in Equ. 1 from Ken nicuttI 
(1998*) then there will be systematic differences between 
the SRFs estimated from the rest-frame UV fiux and the 
true SFRs. 

We have also considered the possibility that the 
Lyman-a emission in one or both of these objects could 
be due all or in part to AGN activity. Given the lack 
of detection of significant emission from highly ionized 
species, such as N[V] and C[IV], we conclude that it is 
unlikely that AGN are playing a significant role in the 
observed Lyman-a emission. The disagreements between 
the two SFR estimates for the two sources discussed here 
highlight the difficulty in making robust SFR measure- 
ments from source- frame UV observables. 

Having accounted for the magnification due to grav- 
itational lensing, we can also compare the lumi- 
nosities of these two sources to large sa mples of 
LAEs and LBGs at co mpara b le redshifts (|Hu et al.l 
20041; iShimasaku et al.' 2006; 'Bouwens et al.' '20071 
Dawson et a l. 2007; Muravama ct al. 2007; Ouchi et al] 
2008t iMcLure et all l2009l ' e.g..). Studies of luminosity 
functions at high redshift are challenging, and parame- 
ters such as L* and a are not nearly so well measured as 
they are in the nearby universe. For comparing Lyman-a 
luminosities, we take the average of 6 different measure- 
ments o^^'Ly-a values from Ouchi et al. (2008) that were 
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fit using samples of LAEs at z ~ 3.7 and z ~ 5.7, assum- 
ing three possible values of the faint end slope of the lumi- 
nosity function, a. For the purpose of comparing contin- 
uum luminosities, we take L^y as the average of the two 
Ltr^nn best-fi t values for LAEs at z ^ 3.7 and z ^ h.l from 
lOuchi et all (HqM)- SGAS J091541-I-382655 is (L2j^_„, 
^Jv) < (0-6L2j,_„, 2Lljy), and SGAS J134331-I-415455 
is (Lly-a; ^hv) = (0-5L2,_„, 0.9L^y). Both of these 
galaxies live at <L* on the luminosity function for sim- 
ilar galaxies at comparable redshifts, which makes them 
interesting targets for studying the properties of typical 
LAEs at z ~ 5 on an individual basis. 

We also attempt to measure the morphologies of both 
lensed LAEs in our Gemini/GMOS imaging and place 
constraints on the intrinsic sizes of these two galaxies. 
Neither source is detected at S/N > 12, which limits 
our ability to make detailed morophological measure- 
ments or shape fits, but we can check the object pro- 
files against the image PSFs and use our estimates of the 
magnification due to lensing to place some rough con- 
straints on the physical sizes of these galaxies. Studies 
of LAE morphologies at z = 3.1 using HST imaging find 
that the Lyman-a and restframe UV emission from these 
sources are spacially coincident an d have half-light radi i 
< l.Skpc (Bo nd et all [2009'. '20H. IBond et al.l (l2009l ) 
argue that S/N> 30 is necessary to make robust mea- 
surements of half-light radii for LAEs, and find that it is 
difficult to distinguish between resolved and unresolved 
compact cores at S/N< 30. The two lensed LAEs are 
best-detected in our Gemini/GMOS i— band imaging, so 
this is the data we use for our morphological measure- 
ments. 

SGAS J134331-H315455 is slightly elongated in the tan- 
gential direction with respect to the center of the cluster. 
We measure its profile to be consistent with a gaussian 
along the tangential axis with a FWHM= 1-1", and con- 
sistent with a gaussian of FWHM= 0.57" along the radial 
axis. The FWHM of the PSF measured from stars in the 
image is ~ 0.57", suggesting that the LAE is unresolved 
along one axis and resolved along the other. This kind of 
morphology is natural for a strongly lensed background 
source where the magnification is generally much greater 
along one axis that the other; a source located on the 
tangential caustic will be highly magnified in the tan- 
gential direction with respect to the center of the lensing 
potential, and magnified very little or not at all in the ra- 
dial direction. We conclude that SGAS J134331-h415455 
is unresolved in the radial direction and barely resolved 
in the tangential direction. Assuming that all of the 
magnification factor of m = 12 is applied in the tan- 
gential direction, we convert a FWHM of 1.1" into a 
physical linear size of '~ 0.6kpc, which indicates that 
most of the emission from this galaxy originates from 
a region very compact in size. This constraint must be 
taken with the caveat that the data are lower S/N than 
would be optimal, but the signficant linear magnification 
of these sources due to gravitational lensing does enable 
us to probe very physically interesting scales even with 
ground-based imaging. We also measure the locations 
of SGAS J13433H-415455 to be identical in the Gem- 
ini/GMOS i— and z— band images to the accuracy of our 
astrometry, 0.2", corresponding to physical scales of 

lOOpc. Based on the observer frame equivalent width 




1 10 

Ll„ (10^= ergs ' h„f) 



Fig. 6. — Two empirical cumulative Lyman-a line emission luminos- 
ity functions (uncorrected for completeness) for LAEs at high redshift 
are plotted; a sample of 2 ~ 4.5 LAEs published by Dawson et al. (2007) 
is giv en by the solid histogram and a sample of 2 ~ 5.7 LAEs (Hu et all 
I2004f) by the dot-dashed histogram. Both the apparent (magnified) and 
intrinsic Lyman-a luminosities for the two sources presented in this 
paper. Open diamonds indicate the apparent luminosities and filled 
squares indicate the intrinsic luminosities, corrected for the magnifica- 
tion (where we have only an upper limit to the intrinsic luminosity for 
SGAS 091541-1-382655). The magnification due to strong lensing allows 
us to probe deep down the Lyman-a luminosity function with sources 
that are observationally the brightest two LAEs known at 2 ~ 5. 

of this source, the fiux in the band filter is ^ 60% 
Lyman-a line emission, whereas the z— band filter mea- 
sures only rest-frame UV continuum emission. Given the 
spatial coincidence of the LAE in the i— and z— bands we 
conclude that the UV continuum and Lyman-a emission 
cannot originate from regions separated by more than 
of order ~ lOOpc projected onto the sky. We have no 
measurement of the relative velocities of the UV con- 
tinuum and the Lyman-a emission and therefore cannot 
constrain a possible separation in velocity. 

SGAS J091541-H382655 is similar to SGAS 
J134331-I-415455 in that the LAE is elongated on 
the tangential axis relative to the center of the lensing 
cluster. Using the GALFIT subtracted i— band image, 
we measure the LAE to be resolved in the tangential 
direction, having a shape consistent with a gaussian of 
FWHM= 1.82", which equates to a constraint on the 
physical scale of < 1.2kpc, where the lower limit on 
the magnification of this source translates to an upper 
limit in the physical size (and again assuming that the 
magnification is entirely along the tangential diretion). 
We also measure the LAE to have a gaussian shape 
with FWHM = 0.75" in the radial direction, which - 
given the FWHM- 0.73" measured from the PSFs of 
stars in the image - is consistent with the source being 
unresolved in the radial direction. In the observer frame 
the equivalent width of the Lyman-a emission for this 
source is 'EWLy-a = 157A, from which we estimate that 
Lyman-a photons are contributing less than 20% of 
the i— band flux, so that we cannot make any attempt 
to measure a spatial offset between the Lyman-a and 
UV continuum emission. Our size constraints for both 
sources are in good agreement with the literature, in 
that we find the sizes of the LAEs to be consistent with 
compact sources of size < l.Skpc. 
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Fig. 7. — The best-fit spectral energy distribution for SGAS 
J091541-t-3S2655 is plotted with a solid line on top of the photometric 
data. The dotted and dashed lines are dust-free, single burst mod- 
els from CB07, with ages of 10 Myr and 50 Myr respectively, scaled to 
match the observed ^-band flux. This illustrates that the photometry 
favors a young (i.e. blue) stellar population, the 50 Myr model already 
deviates considerably from the best-fit SED. Filter transmission curves 
corresponding to the r, i, z and J photometry are plotted in the bottom 
panel. 



1 xlO"'' 



i 4x10"" - 




Fig. 8. — Dust-free single burst CB07 models are scaled to match 
the 2-band flux and plotted on top of the photometric data for SGAS 
J134331-t-415455. Models with ages anywhere in the range of 100 Myr 
to 400 Myr are potentially good fits for the limited data in-hand. Filter 
transmission curves corresponding to the r, i, z and J photometry are 
plotted in the bottom panel. 



3.2. Stellar Mass and UV Continuum Properties 

With several photmetric measurements in hand for 
each of our sources it becomes interesting to investigate 
the spectral energy distributions (SEDs) of the LAEs to 
try and recover the properties (e.g. stellar mass, age, 
SFR, and dust content) of the underlying stellar popu- 
lations in a way that does not rely on as many prob- 
lematic assumptions as, for example, the SFR estimated 
from the rest-frame UV continuum flux that we calcu- 
late in Section 3.1 above. We compare our photometry 
for these sources against modified stellar po pulation syn- 
thesis models. More details can be found in lWuvts et al.l 
(|2OT0f) . we only summarize the main procedure here. We 
use the revised templates of Bruzual and Chariot (CB07, 



based on lBruzual fc CharlotI (|2003f )) with sola r metallic- 
ity, a Chabrier initial mass function (jChabrieii 2003) and 
iCalzettil (|2OOO0 dust extinction law. We investigate a 
range of exponentially declining star formation histories 
of the form SFR{t) ^ exp(t/T), with e— folding times 
r =0.01, 0.05, 0.1, 0.2, 0.5, 1, 2 and 5 Gyr, as well 
as single bursts (SSP) and continuous star formation 
models (CSF). A n updated version of the code Hyperz 
(|Bolzonellall2000D is used to obtain the best-fit SED at 
the fixed spectroscopic redshift of the source via a max- 
imum likelihood procedure. 

The best-fit SED for SGAS J09154H-382655 is shown 
in Figure [7] and corresponds to a dust-free single burst 
with an age of 1.4 Myr. It is important to supplement 
the best-fit stellar population parameters with confidence 
intervals allowed by the photometric uncertainties. We 
create 1000 fake realizations of the observed SED by 
perturbing each broadband magnitude measurement in 
a manner consistent with its errorbars. This set of fake 
SEDs is fit in exactly the same manner as described above 
for the observed SED; bad fits with > 3 are excluded. 
This procedure results in a very young (age < 5 Myr) and 
dust-free {E{B -V) < 0.01) stellar population. Though 
we can place constraints on the age of the observed stel- 
lar population, are data are not sufficient for the SED 
fits to simultaneously constrain the star formation his- 
tory, and consequently the current star formation rate, 
for this galaxy. After correcting for a lensing magnifica- 
tion factor of > 10 , we find a Chabrier stellar mass of 

SGAS J134331+415455 is detected in only 2 bands, 
which is insufhcient to produce a robust best-fit stellar 
population. Instead, we compare the photometry to sin- 
gle burst, dust-free CB07 models of different ages to ob- 
tain some constraints on the age and stellar mass. The 
results are shown in Figure [5] The age of the observed 
stellar population lies in the range between 100 Myr and 
400 Myr, which makes SGAS J134331+415455 signifi- 
cantly older than the very young stellar population seen 
in SGAS J091541+382655. After correcting for a lens- 
ing magnification factor of 12, we constrain the corre- 
sponding stellar masses to be within the range 2 x 10® 
Mq/iq 7 < Mstars < 6 X 10^ MqHqj. All of our photomet- 
ric data for both of these sources sample light blueward 
of the 4000A break in the source frame so that the stel- 
lar population parameters obtained from the SED fitting 
procedures correspond to the most recent episode of star 
formation in each LAE; we have no power to constrain 
stellar mass that may be present in an underlying popu- 
lation of older stars. 

In addition to SED modeling of the photometry of both 
LAEs, we also examine the modest continuum signal that 
is detected for SGAS J091541+382655. Cross correla- 
tion of our CMOS spectrum for this source over a wave- 
length range of 1230-1350A (s ource frame) against the 
composite LBG spectrum from lShaplev et ahl ( 2003[) re- 
sults in a peak at z 5.2, in agreement with the redshift 
measured from the Lyman-a emission. We explicitly ex- 
clude the Lyman-a emission portion of the spectrum in 
this analysis so that the cross-correlation signal that we 
detect is entirely due to low-significance continuum fea- 
tures redward of Ly man-g. The composite LBG from 
iShaplev et all (|2003f) is plotted with our CMOS spectra 
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for both targets in Figures [5] and |31 and visual compar- 
ison of the continuum against the LBG composite spec- 
trum suggests the presence of several strong UV metal 
absorption lines that are observed in well -studied LBGs 
at lower redshifts (e.g. iPettini et aLl l2000'). The features 
are at too low a significance in our data to claim a robust 
detection of individual absorption lines, but the cross 
correlation signal alone is encouraging given our low res- 
olution spectra and limited integration time. Both of 
the sources presented here are excellent candidates for a 
more aggressive spectroscopic followup effort to explore 
the gas properties and stellar metallicity in representa- 
tive low-mass starforming galaxies at z ~ 5. 



star formation. Our analysis of these two LAEs corrobo- 
rates our current understanding of the nature of Lyman- 
a emitting galaxies at high redshift, and the large mag- 
nification of these sources due to gravitational lensing 
makes them excellent candidates for studying the indi- 
vidual properties of galaxies on the faint end of the h^y-a 
and Lf/y luminosity functions at z ^ 5. We encourage 
efforts to followup these sources aggressively on 8- 10m 
class telescopes in order to better study the properties of 
the underlying stellar populations via continuum light. 
These sources are also excellent targets for space-based 
observations, both with current HST instruments and 
with JWST in the future. 



4. SUMMARY AND CONCLUSIONS 

We have identified two lensed Lyman-a emitting galax- 
ies at z ~ 5 near the cores of strong lensing selected 
galaxy clusters. These sources are among the brightest 
galaxies identified at such high redshift, but their intrin- 
sic luminosities are much lower than the observed flux 
due to magnification by the gravitational potential of 
foreground galaxy clusters. We use the available data to 
investigate the underlying stellar populations for these 
galaxies and find that the light - continuum and line 
emission - for SGAS 091541+382655 likely originates 
from a population of young stars with low dust con- 
tent. Both sources are in the process of undergoing active 
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